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Abstract 

We present results from an analysis of B° — ► p + p~ using 316 fb^ 1 of T(45) — ► BB decays observed 
with the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. We measure 
the B° — > p + p~ branching fraction, longitudinal polarization fraction and the CP- violating 
parameters 5i ong and Ci ong : 

B{B°^p + p-) = (23.5 ± 2.2(stat) ± 4.1(syst)) x 10" 6 , 

h = 0.977 ±0.024(stat)+°;^(syst), 

S\on g = -0.19±0.21(stat)+°;^(syst), 

Cong = -0.07 ± 0.15(stat) ±0.06(syst). 

Using an isospin analysis of B — > pp decays we determine the angle a of the unitarity triangle. One 
of the two solutions, a = [74, 117]° at 68% CL, is compatible with the standard model. All results 
presented here are preliminary. 
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1 INTRODUCTION 



In the standard model, charge conjugation-parity (CP) violating effects in the l?-meson system arise 
from a single phase in the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix [1]. Interfer- 
ence between the direct decay and decay after B°B° mixing in B° — ► p + p~ , p^vr^, ir + ir~ results in 
a time-dependent decay-rate asymmetry that is sensitive to the angle a = arg [— V td V t l/V ud V* b ] in 
the unitarity triangle of the CKM matrix. These decays mainly proceed through a 6 ^ uud tree 
diagram. The presence of penguin loop contributions introduces additional phases that shift the 
experimentally measurable parameter a e g away from the value of a. Figure 1 shows the leading or- 
der tree and gluonic penguin loop contributions to this decay. Measurements of the B + — ► p + p° [2] 
and B° — ► p°p° [3] branching fractions show that the penguin contribution in B —>■ pp is smaller 
than the leading order tree diagram [3]. Results from B —> tttt decays [4, 5] and B° —> p^ir 11 [6, 7] 
are discussed elsewhere [8] . 



W+ 




d d d 



Figure 1: Tree and gluonic penguin diagrams contributing to the process B — ► p + p~ . The penguin 
contribution coming from the diagram with a top quark in the loop dominates as contributions from 
processes with u and c quarks are suppressed. 



In B — ► pp decays, a spin particle decays into two spin 1 particles (as shown in Fig. 2). 
Subsequently each p meson decays into a inr pair: p^ — > 7r ± 7r° and p° — > tt + tt^. As a result, 
the CP analysis of B° — > p + p~ is complicated by the presence of one mode with longitudinal 
polarization and two modes with transverse polarization. The longitudinal mode is CP even, while 
the transverse modes contain CP-even and CP-odd states. The decay is observed to be dominated 
by the longitudinal polarization [9], with a fraction fi defined as the fraction of the helicity zero 
state in the decay. Integrating over the angle between the p decay planes (ft, the angular decay rate 
is 



d 2 T 9 



Td cos Bid cos 62 4 



f L cos 2 9 1 cos 2 2 + 1(1 - f L ) sin 2 0! sin 2 2 



(1) 



where #1=1,2 are the angles between the 7r° momentum and the direction opposite to that of the B° 
in the p rest frame. 

In this paper we present an update of previous BABAR measurements of the B° — > p + p~ 
branching fraction, longitudinal polarization fraction, CP-violating parameters and measurement 
of the CKM angle a, previously reported in Ref. [9]. There are several theoretical predictions of the 
branching fraction and /l in B° — > p + p~ decays [10]. The branching fraction of B° — > p + p~ can 
be used to provide constraints on calculations of the branching fractions and charge asymmetries 
in B — > 7T7T and Kit decays [11]. 
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Figure 2: A schematic of the decay of a B meson via two vector particles, p + and p~ , into a four 
pion final state. The p meson final states are shown in their rest frames, and eft is the angle between 
the decay planes of the p mesons. 



The analysis reported here uses 316 fb of data, which is significantly more than our previous 
branching fraction result (82 fb' 1 ) or fi and CP-violating parameter results (210 fb' 1 ). We have 
changed the selection requirements in order to reduce background with only a modest reduction 
in signal efficiency. The new improved probability density-function (PDF) models better account 
for the correlations between variables entering the maximum likelihood (ML) fit, and result in a 
reduced systematic uncertainty of the final results. 

2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric- 
energy B Factory at SLAC. This represents a total integrated luminosity of 316 fb' 1 taken at the 
T(4S) resonance (on-peak), corresponding to a sample of (347 ± 3.9) x 10 6 BB pairs. An additional 
27.2 fb' 1 of data, collected at approximately 40 MeV below the T(4S) resonance (off-peak), were 
used to study background from non-resonant (continuum) e + e~ — > qq events, where q = u,d,s,c. 

The detector is described in detail elsewhere [12]. Surrounding the interaction point is a five 
double-sided layer silicon vertex tracker (SVT) which measures the impact parameters of charged 
particle tracks in both the plane transverse to, and along the beam direction. A 40 layer drift 
chamber (DCH) surrounds the SVT and provides measurements of the momenta for charged par- 
ticles. Both the SVT and the DCH operate in the magnetic field of a 1.5 T solenoid. Charged 
hadron identification is achieved through measurements of particle energy-loss in the tracking sys- 
tem and the Cherenkov angle obtained from a detector of internally reflected Cherenkov light. A 
CsI(Tl) electromagnetic calorimeter provides photon detection, electron identification, and it re- 
construction. Finally, the instrumented flux return of the magnet allows discrimination of muons 
from pions. We use the GEANT4 [13] software to simulate interactions of particles traversing the 
BABAR detector. 
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3 EVENT SELECTION 



We reconstruct B° — > p + p~ candidates (B rec ) from combinations of two charged tracks and two ir° 
candidates. We require that both tracks have particle identification information inconsistent with 
the electron, kaon, and proton hypotheses. The 7r° candidates are formed from pairs of photons, 
each of which has a measured energy greater than 50 MeV. The reconstructed ir° mass must satisfy 
0.10 < m 77 < 0.16 GeV/c 2 . The mass of the p candidates must satisfy 0.5 < m^^o < 1.0 GeV/c . 

Continuum events are the dominant background which is reduced by requiring that | cos(TB, TR)\ 
be less than 0.8, where | cos(TB, TR)\ is the absolute value of the cosine of the angle between the 
B rec thrust axis and that of the rest of the event (ROE) calculated in the CM frame. To distinguish 
signal from continuum we use a neural network (jV) combining the following eight discriminating 
variables. 

• The coefficients, Lq,L2, where these are split into sums over the ROE for neutral and charged 
particles. The coefficients are defined as: 

Lo = (2) 

ROE 

L 2 = ^ Pj |cos(^)| 2 (3) 
ROE 

where pj is the particle momentum and ipj is the angle of the particle direction relative to 
the thrust axis of the B candidate in the center-of-mass (CM) frame. 

• | cos(B, Z)\, the absolute value of the cosine of the angle between the direction of the B and 
z axis in the CM frame, where the z axis is along the electron beam direction. 

• | cos(TB, TR)\ as previously defined. 

• | cos(TB, Z)\, the absolute value of the cosine of the angle between the B thrust and the z 
axis. 

• the sum of the transverse momentum p t of the ROE relative to the z axis. 

Continuum backgrounds dominate near |cos#j| = 1, and backgrounds from B decays tend to 
concentrate at negative values of cos 0j. We reduce these backgrounds with the requirement —0.90 < 
cos6»i < 0.98. 

Signal events are identified kinematically using two variables, the difference AE between the 
CM energy of the B candidate Eb and y/s/2 

AE = E B - (4) 

and the beam-energy-substituted mass 

mES = ^J(s/2 + Pl ■ PB ) 2 /E? - p|, (5) 

where y/s is the total CM energy. The B momentum p# and four-momentum of the initial state 
(Ei,pi) are measured in the laboratory frame. We accept candidates that satisfy 5.25 < tubs < 
5.29 GeV/c 2 and —0.12 < AE < 0.15 GeV. An asymmetric AE selection is used in order to reduce 
background from higher-multiplicity B decays. 
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When multiple B candidates are formed, we select the one that minimizes the sum of (m 77 — 
m n o) 2 where m T o is the ir° mass reported in [14]. In 0.3% the same 7r° mesons are used by multiple 
B candiadtes. We randomly select the candidate entering the fit for such events. 

In order to study the time-dependent decay-rate asymmetry one needs to measure the proper- 
time difference At between the two B decays in the event, and to determine the flavor of the other 
B meson (B tag ). We calculate At from the measured separation Az between the B ICC and B tag 
decay vertices [15]. We determine the B iec vertex from the two charged-pion tracks in its decay. 
The St a g decay vertex is obtained by fitting the other tracks in the event, with constraints from 
the -B rec momentum and the beam-spot location. The RMS resolution on At is 1.1 ps. We only 
use events that satisfy |Ai| < 15 ps and have an error on At of less than 2.5 ps. The flavor of the 
S tag meson is determined with a multivariate technique [16]. The performance of this algorithm 
is summarised in Table 1 for the seven mutually-exclusive tag categories. Events are categorised 
according to decreasing signal purity and increasing mistag probability uj. The categories assigned 
correspond to events with leptons, kaons and pions in the decay products of B tSLg . The Untagged 
category of events are dominated by continuum background, have a mistag probability of 50%, and 
are not used in this analysis. 

Table 1: Tagging efficiency e c , average mistag fraction uj c , and mistag fraction difference Alv c 
between B° and B° tagged events for B° — > p + p~ events in each tagging category. 



Category (c) 


e c 




Auj c 


Lepton 


0.080 ± 0.001 


0.032 ± 0.004 


0.002 ± 0.008 


Kaon 1 


0.114 ± 0.001 


0.053 ± 0.005 


-0.012 ± 0.009 


Kaon 2 


0.176 ± 0.002 


0.152 ± 0.005 


-0.015 ± 0.008 


Kaon-Pion 


0.135 ± 0.002 


0.233 ± 0.006 


-0.018 ± 0.010 


Pion 


0.137 ± 0.002 


0.327 ± 0.007 


0.059 ± 0.010 


Other 


0.095 ± 0.001 


0.412 ± 0.008 


0.042 ± 0.001 


Untagged 


0.266 ± 0.002 


0.500 ± 0.000 


0.000 ± 0.000 



Mis-reconstructed signal candidates or self-cross-feed (SCF) signal may pass the selection re- 
quirements even if one or more of the pions assigned to the p + p~ state belongs to the other B in 
the event. These SCF candidates constitute 50.7% (27.9%) of the accepted longitudinally (trans- 
versely) polarised signal. The majority of SCF events have both charged pions from the p + p~ final 
state, and unbiased CP information. These correct track SCF events are denoted by CT SCF. 
There is a SCF component (13.8% of the signal) where at least one track in B rec is from the ROE. 
These wrong track (WT) events have biased CP information, and are treated separately for the CP 
result. The PDF describing WT events is used only to determine the signal yield and polarization. 
A systematic error is assigned to the CP results from this type of signal event. The total selection 
efficiency for longitudinally (transversely) polarised signal is 7.7% (10.5%). 

4 MAXIMUM LIKELIHOOD FIT 

We obtain a sample of 33902 events that enter an unbinned extended ML fit. These events are 
dominated by backgrounds from qq (81.4%) and BB (16.6%) events. The remaining 2% of events 
are signal. We distinguish between the following components in the fit: (i) correctly reconstructed 
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signal; (ii) SCF signal, split into CT and WT parts; (iii) charm background (6 — > c); (iv) 
charm B° background (b —> c); (v) charmless B° backgrounds; (vi) charmless B ± backgrounds; 
and (vii) continuum background. The dominant B backgrounds come from categories (iii) and 
(iv). The dominant charmless backgrounds are summarised in Table 2. The branching fractions of 
B^ — > (din) 1 * 1 decays have been estimated using the measurements of B° — > afn T [17]. 

Table 2: Dominant charmless backgrounds with assumed or measured branching fraction and the 
number of selected events in the data sample. 



Decay 


Branching fraction (10 <i ) 


Number of events 


afvrT 


39.7 ±3.7 [17] 


94 ±9 


p°p + (long) 


19.1 ±3.5 [2] 


70 ± 13 


a+vr 


20 ±20 


55 ±55 


a?vr+ 


20 ± 20 


45 ±45 


af p T (long) 


24 ±2.5 [18] 


40 ±4 




30 ± 30 [19] 


40 ±40 



We allow the charm and the non-resonant B° — ► p + tt~tt° yields to vary in the fit to data. 
The remaining background yields are fixed to their expected values determined from measured 
branching fractions where available [20]. 

The likelihood function incorporates the following previously defined eight variables to distin- 
guish signal from background: tues, AE, At, J\f, and the m^i^o, and cos#j values of the two p 
mesons. For each of the aforementioned components we construct a PDF that is the product of 
one-dimensional PDFs for each of the variables. The PDFs for all of the components are combined 
to give the likelihood function used in the fit. Table 3 summarises the PDF shapes used for different 
components of the signal. 

The B background AE distributions are described by 3 rd order polynomials, except for non- 
resonant B° — > p + 7r _ 7r° which uses a non-parametric (NP) PDF. The m p distributions for true p 
mesons are parametersised using a relativistic Breit-Wigner, and the fake m p (combinatorial vr^vr ) 
distribution is described using 3 rd order polynomials. The vies distribution of b — ► c backgrounds 
is described using a phase-space- motivated distribution [21]. All remaining background shapes are 
described using NP PDFs. 

The continuum distribution for vies is described by a phase-space- motivated distribution [21]. 
The AE and M shapes are modeled with 3 rd and 4 th order polynomials, respectively. The pa- 
rameters of the tues-, A-E and shapes are allowed to vary in the fit to data. The continuum 
m n ± w o distribution is described using a relativistic Breit-Wigner and a 3 rd order polynomial PDF 
for true and fake p contributions, respectively. The cos B% distribution is described by a 3 rd order 
polynomial. The parameters of the m v ± n a and cos6>j distributions are obtained from a fit to the 
off-peak data. 

The signal decay-rate distribution for both polarizations /+(/_) for B tSLg = B° (B°) is given by 

e -|At|/T 

f±(At) = ±Ssm(Am d At) T C cos(Am d At)] (6) 

where r is the mean B° lifetime, Am d is the B°B° mixing frequency, and S = S\ ong or St ra n and C 
= Ci on g or C tran are the CP-asymmetry parameters for the longitudinally and transversely polarized 
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Table 3: PDFs used to parameterize signal distributions. The abbreviations used are as follows: 
long— longitudinal signal Monte Carlo (MC) simulated data; tran— transverse signal MC simulated 
data; true— true signal (or in the case of background, true p meson); SCF— SCF signal; pN — 
polynomial of order N; G — a Gaussian distribution; CB — a Gaussian with an exponential tail 
which takes the form of [22]; NP — a non-parametric PDF defined as smoothed histogram of MC 
simulated data; RBW — a relativistic Breit- Wigner, SigDT — a distribution of the form of Eq. 6, 
convoluted with a triple Gaussian resolution function. Where indicated, the PDF denoted by pN 
acc. refers to the physical distribution multiplied by a polynomial acceptance function of order N. 



Component 


m E s 


AE 


TV 


At 


cos 


% 


rrip 














Fruc 


Fake 


Fruc 


Fake 


true (long) 


CB 


CB+G 


NP 


SigDT 


p4 acc. 


N.A. 


RBW 


N.A. 


true (tran) 


CB+G 


CB+G 


NP 


SigDF 


p4 acc. 


N.A. 


RBW 


N.A. 


CT SCF(long) (TT) 


CB 


P 2+G 


NP 


SigDF 


p6 acc. 


NP 


RBW 


p3 


CT SCF(long) (TF/FT) 


CB+G 


P 2 


NP 


SigDT 


p6 acc. 


NP 


RBW 


p3 


CT SCF(long) (FF) 


CB+G 


pl 


NP 


SigDF 


p6 acc. 


NP 


RBW 


p3 


WF SCF(long) (TF/FT) 


CB+G 


P 2 


NP 


SigDT 


p6 acc. 


NP 


RBW 


p3 


WF SCF(long) (FF) 


CB+G 


pl 


NP 


SigDT 


p6 acc. 


NP 


RBW 


p3 


SCF(tran) 


CB+G 


GG 


NP 


SigDF 


NP 


NP 


NP 


NP 



signal. The parameters S and C describe B-mixing induced and direct CP violation, respectively. 
S and C for the longitudinally polarized WT signal are set to zero. The At PDF takes into 
account incorrect tags and is convolved with the resolution function described below. Since /l is 
approximately 1, the fit has no sensitivity to either S tTa ,n or Ct ran . We set these parameters to zero 
and we vary them in the evaluation of systematic uncertainties. 

The signal At resolution function consists of three Gaussians (~ 85% core, ~ 14% tail, ~ 1% 
outliers), and takes into account the per-event error on At from the vertex fit. The resolution 
function parameters are obtained from a large sample of fully reconstructed hadronic B decays [15]. 
For WT SCF we replace the 5-meson lifetime by an effective lifetime obtained from Monte Carlo 
(MC) simulation to account for the difference in the resolution. The nominal At distribution for 
the B backgrounds is a NP PDF representation of the MC samples; in the study of systematic errors 
we replace this model with the one used for signal. The resolution for continuum background is 
described by the sum of three Gaussian distributions whose parameters are determined from data. 

5 RESULTS 

From the extended maximum likelihood fit described above, we obtain the following results 

615±57(stat), 
0.977 ±0.024(stat), 
-0.19±0.21(stat), 
-0.07±0.15(stat), 

after correction for a +28 event fit bias on the signal yield and a correction for a —0.007 fit bias 
on fi. We discuss the origin of these fit biases in Section 6. The B° — > p^tt^tt background 



N (signal) = 

h = 

S'long = 

Clong = 
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yield obtained from the fit is 9.2 ± 53.6 events. Figure 3 shows distributions of rriES, AE, cos#j 
and m n ± n o for the highest purity tagged events with a loose requirement on J\f. The plot of rriES 
contains 15.6% of the signal and 1.1% of the background. For the other plots there is an added 
constraint that uies > 5.27GeV/c 2 ; these requirements retain 13.9% of the signal and 0.4% of the 
background. Figure 4 shows the At distribution for B° and B° tagged events. The time-dependent 
decay-rate asymmetry 

N(At) - N(At) 



anp(At) _ 
V ' N(At) + N(At) 

is also shown, where N (N) is the decay-rate for B° (B°) tagged events. 
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Figure 3: The distributions for the highest purity tagged events for the variables uies (a), AE (b), 
cosine of the p helicity angle (c) and m^i^-o (d). The dashed lines are the sum of backgrounds, and 
the solid lines are the full PDF. 



6 SYSTEMATIC ERROR STUDIES 



Table 4 lists the possible sources of systematic uncertainties on the values of the signal yield, fi, 
S\ ong and Ciong that have been studied. The systematic uncertainties, listed in Table 4, are briefly 
discussed in the following. 

• The uncertainty from PDF parameterisation is obtained by varying PDF shape parameters by 
±lcr, in turn. The deviations obtained are added in quadrature to give the quoted uncertainty 
from the PDF parameterisation. 

• The systematic uncertainty on the yield from the fraction of SCF events is obtained taking 
the difference between the nominal result and the number of events fitted (after correcting for 
the fraction of SCF events in MC simulated data) when using the true signal PDFs to extract 
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Figure 4: The At distribution for a sample of events enriched in signal for (a) B° and (b) B tagged 
events. The dotted lines are the sum of backgrounds, the dashed lines are the qq background and 
the solid lines are the sum of signal and backgrounds. The time-dependent CP asymmetry (see 
text) is shown in (c), where the curve is the measured asymmetry. 



the yield. The uncertainty on the other signal observables comes from varying the fraction of 
SCF events by 5% per tt° (10% total), which is twice the observed data/MC difference seen 
in our control sample of B° — > D~ p + events. 

• The uncertainty on the rriEs and AE widths is obtained from the observed shifts relative to 
our nominal result, when allowing these parameters to vary independently in the fit to data. 

• We vary the B background normalisation within expectations for each background in turn. 
The deviations obtained when doing this are added in quadrature to give the quoted uncer- 
tainty from this source. 

• As the branching fractions of some of the B backgrounds are not well known, we assign an 
additional uncertainty coming from the maximum shifts obtained when allowing each of the 
fixed backgrounds to vary in turn in the fit to data. 
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• Additional uncertainties on the CP results come from possible CP violation in the B back- 
ground. We use existing experimental constraints where possible, otherwise we allow for a 
CP asymmetry up to 10% in B decays to final states with charm, and up to 50% in B decays 
to charmless final states. 

• The physics parameters r = 1.530 ± 0.009 ps and Anid = 0.507 ± 0.005 Zips [14] are varied 
within the quoted uncertainty. 

• The tagging and mistag fractions for signal and the B backgrounds are corrected for data/MC 
differences observed in samples of fully reconstructed hadronic B decays. Each of the tagging 
and mistag parameters is varied in turn by the uncertainty from the correction. The deviations 
obtained when doing this are added in quadrature to give the quoted uncertainty from this 
source. 

• Allowing for possible CP violation in the transverse polarization, and in the WT longitudinally 
polarised signal SCF events results in additional uncertainties on signal observables. We vary 
S and C by ±0.5 (±1.0) for the transverse polarisation (WT SCF). 

• Possible CP violation from interference in doubly Cabibbo-suppressed decays (DCSD) on the 
tag side of the event [23] contribute to systematic uncertainties on Si on g and Ci on g • 

• We estimate the systematic error on our results coming from neglecting the interference 
between B° — ► p + p~ and other 47r final states: B — ► a\ir, pirir and B — ► iririr°ir°. Strong 
phases are varied between —180 and 180 degrees, and the CP content of the interfering 
amplitudes are varied between zero and maximum using uniform prior distributions, and the 
RMS deviation of the parameters from nominal is taken as the systematic error. We assume 
that the amplitudes for B — ► pinr and 7r7T7r 7r° are equal to the amplitude of B — > ain when 
calculating this systematic uncertainty. 

• As the PDFs used in the ML fit do not account for all of the correlations between discrim- 
inating variables used in the fit, the results have a small bias. We calculate the fit bias 
on the signal observables from ensembles of experiments obtained from samples of signal 
and charmless B background MC simulated events combined with charm and qq background 
events simulated from the PDFs. The 28 event bias on signal yield and 0.007 bias on fi is 
corrected, and 100% of the calculated fit bias is assigned as a systematic uncertainty on all 
signal observables. We do not observe a significant bias on S'long an d Ci ong . 

• Small imperfections in the knowledge of the geometry of the SVT over time can affect the 
measurement of S'long and Ci on g- We vary the alignment according to the results obtained 
from the study of e + e~ — ► e + e~,p + p~ events in order to estimate the magnitude of this 
systematic error on our CP results. 

The branching fraction has multiplicative systematic uncertainties from the reconstruction of ir° 
mesons in the detector (6%), uncertainties in the reconstruction of charged particles (0.8%), and the 
discrimination of tt^ from other types of charged particles (1%). In addition to these uncertainties, 
there is a 1.1% uncertainty on the number of BB pairs in the data sample. Statistical uncertainties 
arising from the MC samples used in this analysis are negligible. 
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Table 4: Summary of additive systematic uncertainty contributions. 



Contribution 


0"(A ' signal) 




cr (5'long) 


c(Clong) 


Pub parameterisation 


+16.7 
-30.2 


+0.0082 
-0.0064 


+0.0149 
-0.0425 


+0.0300 
-0.0306 


SCF fraction 


84.0 


+0 0007 
-0.0011 


-0.00355 


_i_n 0070 
-0.00683 


m,ES and AE width 


22.9 


0.005 


0.011 


0.012 


B background normalisation 


-1-1 fi 

-17.2 


_i_n nn^^ 
-0.0038 


_i_n nno.fi 
-0.0115 


-0.0015 


tI nQ "f" i n it h\ q p Lrmr^i nine 


oo.u 




u.uoo 


u.uuu 


CPV in B background 


+3.3 
-2.0 


+0.0006 
-0.0016 


+0.0059 
-0.0214 


+0.0118 
-0.0115 


r 


+0.1 
-0.4 


i r\ r\r\r\r\ 

+0.0000 
-0.0002 


+0.0002 
-0.0008 


0.0007 


Am 


+0.0 


+0.0000 


+0.0014 


+0.0018 


-0.2 


-0.0002 


-0.0020 


-0.0012 


tagging and dilution 


+2.6 
-8.1 


+0.0029 
-0.0021 


+0.0016 
-0.0053 


+0.0068 
-0.0054 


transverse polarisation CPV 


+0.0 

-8.3 


+0.0057 
-0.0000 


+0.0125 
-0.0152 


+0.0095 
-0.0110 


WT SCF CPV 


+0.2 
-1.1 


+0.0000 
-0.0003 


+0.0051 
-0.0065 


+0.0116 
-0.0113 


DCSD decays 






0.012 


0.037 


Interference 


14.8 


0.0036 


0.023 


0.022 


Fit Bias 


28 


0.007 


0.002 


0.022 


SVT Alignment 






0.0100 


0.0055 


Total 


+97 
-101 


+0.015 
-0.013 


+0.05 
-0.07 


±0.06 
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7 FINAL RESULTS 



Our results are 



B(B°^p + p-) = (23.5 ± 2.2(stat) ± 4.1(syst)) x 1(T 6 , 




where the correlation between Si ong and Ci ong is -0.058. This measurement corresponds to a e s = 



and CP parameters are in agreement with our earlier publication [9], with significantly improved 
precision. 

We constrain the CKM angle a from an isospin analysis [24] of B — ► pp. The inputs to the 
isospin analysis are the amplitudes of the CP-even longitudinal polarization of the pp final state, 
as well as the measured values of Si on g and Ci ong for B° — > p + p~. We use the following numerical 
inputs in the isospin analysis: 

• The average of the measurements of B(B° — ► p + p~) and /l, presented here, and those 
reported in Ref. [25]. 

• The combined branching fraction and Jl for B — > p + p° from Ref. [2]. 

• The central value corresponding to the upper limit of B(B — > p°p°) from Ref. [3]. 

• Siong and Ci on g presented here. 

We ignore possible 1 = 1 amplitudes [26] and electroweak penguins in this paper. 

To interpret our results in terms of a constraint on a from the isospin relations, we construct a 
X 2 that includes the measured quantities expressed as the lengths of the sides of the isospin triangles 
and we determine the minimum Xq- We have adopted a simulated experiment technique to compute 
the confidence level (CL) on a; our method is similar to the approach proposed in Ref. [27]. For 
each value of a, scanned between and 180°, we determine the difference ^Xdata( q; ) between the 
minimum of x 2 ( a ) an d Xo- We then generate MC experiments around the central values obtained 
from the fit to data with the given value of a and we apply the same procedure. The fraction of 
these experiments in which Axjy[ C (a) is smaller than Axp ATA (a) is interpreted as the CL on a. 
Figure 5 shows 1 — CL for a obtained from this method. Selecting the solution closest to the CKM 
combined fit average [28, 29] we find a = [74, 117]° at 68% CL, where the error is dominated by 
\ a e& ~ a \ which is 18° at 68% CL. The constraint obtained on a is less precise than the previous 
result in Ref. [9] for the following reasons: (i) the improved world average branching fraction for 
B+ -> p+p° has come down, and (ii) the results of the latest search for B° — > p°p° show evidence 
for a signal. Both of these factors lead to an increase in the the penguin contribution to the total 
uncertainty on a when using an isospin analysis. 

It is possible to calculate a model-dependent constraint on a using SU(3) as detailed in Ref. [30]. 
This model relates the penguin amplitude in B + — > K*°p + to the penguin amplitude in B° — > p + p~ , 
and can be used to obtain a more precise, but model dependent constraint on a. 




The measured branching fraction, polarization, 
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Figure 5: Confidence level on a obtained from the isospin analysis with the statistical method 
described in [28]. The dashed lines correspond to the 68% (top) and 90% (bottom) CL intervals. 



8 SUMMARY 

In summary we have improved the measurement of the branching fraction, fi, and CP-violating 
parameters Si ong and Ci ong in B° — ► p + p~ using a data-sample of (347±3.9) x 10 6 BB pairs. We do 
not observe mixing-induced or direct CP violation. We derive a model-independent measurement 
of the CKM angle a. The results presented in this paper are consistent with previous results 
presented by BABAR [9] and Belle [25]. 
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